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A distance estimate based on angular expansion 
for the planetary nebula NGC 6881 
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ABSTRACT 

In this paper, we report on high angular resolution radio observations of the planetary nebula 
NGC 6881 obtained with the Very Large Array at a wavelength of 6 cm. The emission appears 
to be the superposition of a roundish core and a point-symmetric bipolar structure elongated 
along a position angle of about 145°. This is strongly reminiscent of the morphology seen in 
Ha and [Nil] images. A comparison between VLA observations obtained in 1984 and 1994 
clearly reveals the expansion of the core of the nebula, at a rate of 2. 1 + 0.7 mas yr~' . Assum- 
ing that the expansion velocity in the plane of the sky (determined from these measurements) 
and the expansion velocity along the line of sight (estimated from optical spectroscopy avail- 
able in the literature) are equal, we find a distance to NGC 6881 of 1.6 ± 0.5 kpc + 0.3 kpc, 
where the first error reflects the uncertainty on the expansion, and the second error comes 
from the potential difference between pattern and material speeds. This distance is compat- 
ible with (but does not necessarily imply) an association of NGC 6881 with the nearby HII 
region Sh 2-109 and, more generally, the Cygnus star-forming region. 

Key words: planetary nebulae: general — planetary nebulae: individual (NGC 6881, 
PG 74.5+02. 1) — astrometry — stars: late type 



1 INTRODUCTION 

NGC 6881 (PG 74.5+02.1, IRAS 20090+3715) is a very in- 
teresting planetary nebula (PN) with a system of multiple 
bipolar lobes trace d by the emission from ionized as well 
as neutral material ( Guerrero & Manchado 19981: Guerrero et"al] 



200d: iKwok & Sul 12005]: 



Ramos-Larios et al 



200i) 



, Originally, 

Guerrero & Manchadol d 19981) classified it as a quadrupolar plan- 



etary nebula based on the examination of the [N II] emission and 
the [N II]/[0 III] and [N Ill/Ha ratio images. These authors found 
one pair of lobes consisting of two highly collimated structures that 
finish as bright knots where the emission of [Nil] is enhanced (outer 
lobes). With a size of about 8", the inner lobes are shorter but much 
brighter than their outer counterparts. Moreover, they exhibit a clear 
point-symmetric morphology. The expansion velocity derived for 
these lobes is 90 km s -1 , which yields a kinematical age of about 
400 rf years (where d is t he source distance in kpc). 

IKwok & Sil J2005h obtained Hubble Space Telescope images 
of this nebula from which they identified one more pair of ion- 
ized lobes. They have the same extension and morphology as the 
inner lobes but with a position angle slightly different. This dis- 
crepancy is not well understood. These authors also found several 
rings circumscribed in the inner walls of the ionized bipolar lobes. 



The presence of these multipolar lobes and rings oriented in dif- 
ferent directions suggests the action of precessing collimated fast 
outflows that change t h eir pr ecession axis as a function of time. 
I Guerrero & Manchadol d 19981) proposed this to explain the exis- 
tence of a loop-like structure traced by the ionized emission, lo- 
cated at the southeast edge of the lobes and elongated in the direc- 
tion perpendicular to the lobes. They modeled the loop as a pre- 
cessing jet that , assuming ballistic expansi on, started ~7400rf years 
ago. However, iRamos-Larios et all d2008l) found that the location 
of the loop-like structure coincides perfectly with the ending tip 
of the southeast H2 lobe, suggesting that it is marking a region of 
interaction of the bipolar lobe with dense material. These authors 
noted that NGC 6881 is projected towards the HII region Sh 2-109 
and suggested that they are physically associated. They proposed 
that the southeast lobe of NGC 6881 is interacting with dense ma- 
terial, forming the observed loop-like structure and the sharp edge 
of the H2 lob e. The HII region S h 2-109 is located at a distance of 
1. 4+0.4 kpc dFich & Blitzll 19841). Using a distanc e scale proposed 
bv lCahn et all dl992h . lRamos-Larios etal1d2008r) calculated a dis- 
tance of ~1.5 kpc to NGC 6881, which supports the asso ciation 
with t he region Sh 2-109. However, the scale developed by IZhand 
d 19951) yields a distance of 6.4 kpc, in which case they would not 
be physically associated. 
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Statistical techniques can provide distance estimates to PNe 
as accurate as 50% on average, but when they are applied to in- 
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Table 1. Observational parameters 



Epoch 


Phase Calibrator 


Bootstrapped Flux density 


1984 Dec 29 (1984.99) 


2023+336 


2.07 ±0.01 Jy 


1994 Mar 31 (1994.24) 


2005+403 


2.98±0.01 Jy 



dividual objects the errors can be as large as a factor of 2 or 
even more (e.g. iTafova et al.l 1201 fl) . This results in a very un- 
certain estimation of the Galactic location and physical parame- 
ters of PNe. In order to obtain more accurate estimations of the 
distance to PNe, other techniques have been used. One that has 
proven effective is the so-called expansion parallax method, which 
is based on the comparison of the angular expansion of the source 
on the plane of the sky, with the velocity expansion along the 
line of sight meas ured from the width of some appropriately cho- 
sen spectral lines dMassonlll986h . In principle, the angular expan- 
sion could be measured using any tracer, but the free-free radio 
emission fr om the ionized gas h as been, by far, the most popu- 
lar choice J"MassorJll986l. 1 19891: iGomez et alj|l993l: IHaiian et all 



19931 119951: iKawamura & Massed 1 199d: IHaiian & Terzianlll996l: 
Christianto & Seaquistl 1 19981 : IGuzman et al.l 120061 : IZiilstra et al.l 
20081 : IGuzman et al.ll20Qg|) 7 Hubble Space Telescope observations 



at optical wave lengths have also occasionally been used to measure 
the expansion dReed et al.lll999l : |Palen et al]2002T) . 

In this work we have used the expansion parallax method to 
measure the distance to the PN NGC 6881. Our new distance esti- 
mate will be used to determine with improved accuracy the physical 
parameters of the source, particularly its size and age. It will also 
enable us to examine the association of NGC 6881 with the HII 
region Sh 2-109. For this purpose, we have retrieved data from the 
Very Large Array (VLA) archive of the National Radio Astronomy 
Observator>Q (§2). The structure of NGC 6881 and the details of 
the component that was used for the measurements are described 
in §3. The technique used to determine the angular expansion is 
presented in §4. The estimation of the distance and the location of 
NGC 6881 in our Galaxy, as well as other physical parameters, are 
discussed in §4.1. 



2 THE DATA 

We have made use of two data sets of observations toward the PN 
NGC 6881 retrieved from the VLA archive. The observations were 
performed in the continuum mode (bandwidth = 50 MHz) at a fre- 
quency of 4.8 GHz (corresponding to a wavelength of 6 cm). The 
configuration of the array for both data sets was the most extended 
one (configuration A). The two data sets correspond to observa- 
tions carried out on 1984, December 29 (1984.99; project AK1 13) 
and 1994, March 31 (1994.24; project AH509), respectively. This 
gives a time baseline of 9.25 years between the observations. The 
data were edited and calibrated using the Astronomical Image Pro- 
cessing System (AIPS) following standard procedures (see Table 
1). Both data sets were self-calibrated in phase and amplitude. 
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Figure 1. Contour image of the 4.8 GHz (6 cm) continuum emission from 
NGC 6881 for the 1994.24 epoch. The contours are -1, 1, 2, 3, 4, 5, 10. 20. 
30, 40, 50, 60, 70, 80, and 90% the peak flux of 9.8 mjy beam" 1 . The rms 
noise of the image is 0.025 mJy beam -1 . The synthesized beam (0"57 X 
0"44 with a position angle of 87°) is shown in the bottom left corner of 
the image. The cross marks the position adopted as the center of the nebular 
emission: o-(2000)= 20* 10'" 52f 446, 5(2000)= +37° 24' 42"43. 



3 THE STRUCTURE OF THE PN NGC 6881 

In the radio continuum images, NGC 6881 clearly exhibits a bipo- 
lar, point- symmetric structure elongated along a position angle of 
145° (see Fig. 1). This morphology is ver y similar to that of the 
Ha, [N il] and [OIII] images presented by I Guerrero & Manchadol 
Jl998l) . although the radio continuum emission image is not as sen- 
sitive in the outermost portions of the lobes as those in Ho- and 
[Nil]. In consequence, neither the knots nor the loop-like structure, 
which are detected in [Nil] light, are detected in radio continuum 
emission. In the central region of the nebula, a two-peaked bright 
core (also seen in [Nil] images) elongated along a position angle 
of 55° is clearly present. Guerrero & Manchado (1998) interpreted 
this structure as an expanding ring in the equatorial plane of the 
nebula. 

If we assume that the expansion velocity of the central core 
along the line of sight is the same as that on the plane of the sky, 
the angular expansion rate should be about 2 mas yr if the source 
is at 1.5 kpc and about 0.5 mas yr~' for d = 6.4 kpc. Such values, 
although small, are measurable, provided the core of the emission 
can be separated from the lobes. For interferometric data, the struc- 
tures associated with different angular scales can be separated by 
restricting the portion of the (u - v) plane considered during the 
imaging process (Gomez et al. 1993). In the case of NGC 6881 the 
emission from the extended lobes can be almost entirely suppressed 
by selecting only the visibilities corresponding to (u - v) spacings 
longer than about 30 kA (corresponding to ~7"). In addition, the ex- 
pansion is best detected if the longest (u—v) spacings are somewhat 
down- weighted and the angular resolution is degraded so the source 
is only a few resolution elements across (e.g. Gomez et al. 1993). In 

1 The NRAO is operated by Associated Universities, Inc. under a coopera- 
tive agreement with the National Science Foundation. 



A distance estimate for the PN NGC 6881 3 





1984.99 



20 10 52.7 52.6 52.5 52.4 52.3 52.2 

RIGHT ASCENSION (J2000) 




1994.24-1984.99 



52.5 52.4 52.3 

RIGHT ASCENSION (J2000) 





1994.24 



52.6 52.5 52.4 52.3 

RIGHT ASCENSION (J2000) 





MODEL 



52.5 52.4 52.3 

RIGHT ASCENSION ( J2000) 



Figure 2. Top: Reconstructed contour images of the 6 cm continuum emission from NGC 6881 for 1984.99 (left) and 1994.24 (right). The contours are -4, 
4, 7, 10, 15, 30, 50, 70, 100, 130, 150, 200, 250, and 280 times 0.07 mJy beam -1 , the average rms noise of the images. The individual rms of the images is 
0.09 mJy beam" 1 for the 1984.99 image and 0.05 mJy beam for the 1994.24 image. Bottom: contour images of the 6 cm difference image (left) and of the 
best model (right) obtained as described in the text. The contours are —11, -9, -7, -5, -3, 3, 5, 7, 9, and 11 times 0.095 mJy beam -1 , the rms noise of the 
difference image. The restoring beam (0"87 X 0"80 with a position angle of 0°) is shown in the bottom left corner of each image. 



the present case, this was achieved by considering only visibilities 
in the (u - v) range from 30 to 550 kA, and by applying a Gaussian 
taper function with a full width at half maximum of 200 kA. Finally, 
the images were reconstructed to minimize possible differences in 
the u-v coverage using a beam of 0!'87 x 0''80, P.A.= 0°, the aver- 
age of the individual beams of each observation. These procedures 
yielded noise levels of ~ 0.09 and 0.05 mJy beam -1 , for the 1984 
and 1994 data sets, respectively. The difference image between the 
two epochs (1994.24-1984.99), was produced following Guzman 
et al. (2006, 2009), and is shown in the bottom left panel of Fig- 
ure 2. This image shows the typical negative and positive pattern 
expected when expansion is present. 



4 ESTIMATION OF THE ANGULAR EXPANSION 

The images of NGC 6881 obtained at epochs 1984.99 and 1994.24 
(upper panels of Figure 2) look very similar, so the expansion of the 
nebula is not immediately obvious. The negative-positive pattern 
seen in the difference image (1994.24-1984. 99), however, clearl y 
reveals the expansi on of the ionized nebula. [Gomez et al I dl993h . 
and more recently iGuzman et al. I d2006l \200% have developed a 



technique that has proven effective to measure very small expan- 
sion rates using observations at two epochs. The technique is based 
on the assumption that the expansion of the ionized nebula is self- 
similar. To characterize the expansion, the image corresponding to 
the better of the two data sets is self-similarly expanded by a factor 
1 + e (if the best data set corresponds to the older observation) or 
shrunk by a factor 1 —e (if it is the more recent), and subtracted from 
itself. This provides a model difference image that can be compared 
with the true difference image. The correct expansion rate clearly 
corresponds to the value of e for which the model difference image 
is most similar to the measured difference image. Quantitatively, 
this corresponds to the value of e that minimizes the y 2 defined as 

X 1 = X , V/ - Tuf/o 2 , (1) 

where Mq is the value at pixel (i, j) in the model difference image, 
Tjj is the value at the corresponding pixel of the true difference 
image and <x is the standard deviation of the data. The sum is taken 
over the appropriate portion of the images. 

In the present case, the data obtained at epoch 1994.24 have 
better quality, so the corresponding image was contracted by a fac- 
tor (1 - e). The value of the x 2 as a function of e is shown in Fig- 
ure 3. It has a clear minimum at e = 0.027. The model difference 
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Figure 3. x 1 fit of the residual image obtained from subtracting the model 
to the data as function of the contraction factor e. The minimum of the fitted 
curve (dashed-line) indicates the value of e that minimizes the^ 2 function. 



Table 2. Physical parameters for NGC 6881 



Parameter 


Value 


Flux (4.8 GHz) [mJy] 


120.0 ± 0.4 


Distance [kpc] 


1.6 ±0.5 


Radius of the central ring [AU] 


1170 ± 380 


Expansion rate of the ring [km s ] 


16 ±2 


Kinematic age of the ring [years] 


350 ± 110 


n e in the central region [10 4 cirT 3 ] 


5 ± 1 


Mi in the central region [10~ 2 M Q ] 


0.22 ± 0.08 


EM in the central region [10 6 cirT 6 pc] 


28 ± 1 



image corresponding to that value of e is shown in the bottom right 
panel of Figure 2. Once the minimum value of e is determined, the 
confidence limits on this value are estimated as follows. The num- 
ber of independent data points in the fit, n, is estimated by dividing 
the solid angle of the region used to estimate e by the solid angle 
of the synthesized beam. Here, the region of the images considered 
to calculate the^ 2 corresponds to n = 20 beams. Assuming that the 
plot in Figure 3 r epresents a normal x" distribution and following 
IPress et all ( 1 1992b . we scale by multiplication the value of the min- 
imum to be 20. The confidence limits are estimated by determining 
the e interval given for a value in the horizontal axis of n + 1 = 21 
d Avn ill 1 976l : [Walll 1 99fj) . This yields e = 0.027 ± 0.010. 

The good agreement between the images produced for the real 
difference and the model difference suggests that the assumption of 
self-similar expansion for the central ring of the PN NGC 6881 
is reasonable, at least in the equatorial plane of the nebula. The 
structure in the direction of the bipolar lobes appears to be more 
complex, and a different expansion velocity (~90 km s~'; Guer- 
rero & Manchado 1989) might be more appropriate. Moreover, as 
we mentioned earlier, the outflow direction may be changing with 
time (Kwok & Su 2005). Therefore, a 2D or 3D model (such as 
that presented by Zijlstra et al. (2008) for NGC 7027) would, in 
principle, be needed. However, constructing a detailed kinematic 



model of NGC 6881 would require unavailable additional observa- 
tions and, at any rate, the value of e found here is almost entirely 
constrained by the expansion of the central core/ring. Thus, as long 
as we restrict ourselves to this specific component, our determina- 
tion provides an accurate estimate of the expansion. From the value 
of e estimated earlier, the angular expansion rate of NGC 6881 can 
be calculated as: 



(2) 



where, 8 is the radius of the ring and At is the time baseline be- 
tween the observations. From the separation of the two peaks seen 
in the full free-free emission images (Figure 1), we estimate that 
the radius of the ring is 0"73 ± 0"01, and an angular expansion rate 
of 6 = 2.1 ± 0.7 mas yr 1 for the central ring of NGC 6881. 



4.1 The distance to NGC 6881 

The distance d to NGC 6881 can be obtained from its angular ex- 
pansion rate as: 

i 



211 



km s 



-i 



e 



mas yr 



(3) 



d 
pc 

where v cxp is the expansion velocity of the central ring along the 
line of sight. As mentioned previously, the expansion velocity of 
the ring has been estimated fro m the line splitting of the [Nil] by 
iGuerrero & Manchadol d 19981) to be v exp = 14.0 km s _1 . Other es- 
timates of the expansion velocity of NGC 6881 using [ QUI] lines 
resulte d in values of 18 km s _I dRobinson et alj[l982l) and 16.5 
km s 1 dWeinbergerll989l) . Thus, we will use an average expansion 
velocity value of v exp = 16 ± 2 km s _1 . This leads to a distance d = 
1.6 ± 0.5 kpc; table 2 shows the physical parameters of NGC 6881 
calculated using this value. 

Distance measurements based on angular expansion can po- 
tentially be affected by the difference between the pattern velocity 
(measured by the expansion) and the true gas velocity (measured 
by the radial velocities). Mellema (2004) discusses this effect in 
detail, and concludes that it typically results in errors on the dis- 
tance of about 20%. This would add a "systematic" error to our 
measurement of 0.3 kpc. We conclude that the distance to NGC 
6881 is 1.6 ± 0.5 ± 0.3 kpc, where the first error reflects the un- 
certainty on the expansion, and the second error comes from the 
potential difference between pattern and material speeds. 

The distance found here is in g ood agreement with that ob- 
tained by iRamos-Larios et aT I d2008h . d ~ 1.5 kpc. It is also com- 
patible with the physical association of NGC 6881 with the nearby 
HII region Sh 2-109 (see §1). When s uperimposed over the ex- 
tended HI emission jCappaetal.il 19961) . NGC 6881 appears pro- 
jected toward a cavity with a vlsr ~ 10 km s" 1 , very close to 
that of the Cyg OB 1 association. The esti mated distance to this as- 
sociation lies in the r ange 1.25-1.83 kpc dLozinskava et al.ll 19981 : 
ISchneider et al.ll2006l) .The association Cyg OB 1 is located to the 
south-west of the Cygnus star-forming region. Recently, Rygl et 
al. (2010) measured very accurately the parallax of the massive 
star-forming region W75N with the Cygnus region and found d = 
1.32^ 'Jo kpc. This suggests that the entire Cygnus complex (which 
is about 200 pc across) is located to distances between 1.2 and 
1.8 pc. Thus, our result supports the idea that NGC 6881 might 
be associated with the massive star forming region in Cygnus. 
Although this may be a simple coincidence, it may also indicate 
that the progenitor of NGC 6881 could have been a short-lived 
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intermediate-mass star. Interestingly, some of the observed char- 
acteristics of NGC 6881 m ay be readily interpreted in this scheme. 
iRamos-Larios et aD d2008h observed that the southeast H2 lobe of 
NGC 6881 is less extended than its northwest counteipart, and 
shows a sharp edge. They note that a possible explanation for this 
asymmetry would be an interaction between the lobes of NGC 6881 
and inhomogeneous interstellar gas. This could occur naturally if 
NGC 6881 were indeed associated with the star-forming region. 



5 CONCLUSIONS 

We have presented observations of the planetary nebulae 
NGC 6881 obtained using the VLA at 6 cm (4.8 GHz) at two 
epochs separated by 9.25 years. Assuming a self-similar expan- 
sion for the ionized gas, we determine an expansion angular rate 
for NGC 6881 of 2.1 ±0.7 mas year 1 , and a distance of 1.6 ± 0.5 
kpc ±0.3 kpc, where the first error reflects the uncertainty on the 
expansion, and the second error comes from the potential differ- 
ence between pattern and material speeds. This places NGC 6881 
at the same distance as the HII region Sh 2-109 and the Cygnus 
star-forming region. 
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